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Abstract 

The SU{2) gauge invariant Dirac- Yang-Mills mechanics of spatially homo- 
geneous isospinor and gauge fields is considered in the framework of the 
generalized Hamiltonian approach. The unconstrained Hamiltonian system 
equivalent to the model is obtained using the gaugeless method of Hamil- 
tonian reduction. The latter includes the Abelianization of the first class 
constraints, putting the second class constraints into the canonical form 
and performing a canonical transformation to a set of adapted coordinates 
such that a subset of the new canonical pairs coincides with the second 
class constraints and part of the new momenta is equal to the Abelian 
constraints. In the adapted basis the pure gauge degrees of freedom au- 
tomatically drop out from the consideration after projection of the model 
onto the constraint shell. Apart from the elimination of these ignorable 
degrees of freedom a further Hamiltonian reduction is achieved due to the 
three dimensional group of rigid symmetry possessed by the system. 
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I. INTRODUCTION 



The correct canonical formulation of the quantum theory of non-Abelian fields assumes 
a detailed knowledge of the corresponding classical generalized Hamiltonian dynamics [0]- 
Since the introduction of non-Abelian gauge fields by C.N. Yang and R.L. Mills [^] over 
forty years ago essential progress in this direction has been made. Rigorous statements 
about the geometrical structure of the configuration and the phase space have been estab- 
lished. It turned out that due to the underlying non-Abelian gauge symmetry the "true 
phase space" of Yang-Mills theory, namely the quotient space of phase space by the action 
of gauge transformations, possesses a rich topological structure In the framework of 
traditional perturbation theory these geometrical peculiarities are not taken into account 
and as a result the description of large scale effects, including confinement, is beyond its 
scope. The most important lesson one has learned is that, in order to reach a complete de- 
scription, it is necessary to first reformulate Yang-Mills theory in terms of gauge-invariant 
variables and only after this step apply any approximation method. With this aim sev- 
eral different representations for the physical degrees of freedom of non-Abelian theories 



T9[| have been proposed. All these approaches lead to an unconstrained Hamiltonian 



system, which exhibits non-perturbative features and are in some sense alternatives to 
the conventional perturbative approach. Whereas perturbation theory is appropriate for 
the computation of short distance effects, the unconstrained formulation is adapted to 
the study of large scale phenomena if the gauge invariant expressions are evaluated in a 
derivative expansion. Since the work by Matinyan et al. pO[, the corresponding zeroth 
order or long-wavelength approximation, the Yang-Mills mechanics of spatially homoge- 



neous gauge fields, has been studied extensively from different points of view (see e.g. [21 



2^ and references therein). In the present note we shall continue the study of the model 



arising in this approximation, pursuing the aim to prepare the necessary background for 
studying the problem of construction of the reduced phase space of QCD. Due to the 
spatial homogeneity condition conventional Dirac- Yang-Mills theory reduces to a theory 
describing a finite dimensional system which is incomparably simpler than the exact field 
system. At the same time, however, it possesses all the main peculiarities of the full 
theory and can be used as a laboratory for testing the viability of ideas and techniques 
that could be applied in the general case. 

Below we shall isolate the true dynamical degrees of freedom of SU (2) Dirac- Yang- 
Mills theory in the long- wavelength approximation using the gaugeless approach to the 
reduction in the number of degrees of freedom instead of the conventional gauge fixing 
method. The cornerstones for this method applied to a system with first class con- 
straints are the procedure of Abelianization of constraints (replacement of the original 



Presumably, S.Shanmugadhasan |25| was the first to employ the classical Lee-Cartan method 
of reduction (see e.g. [26| - |^]) in the framework of generalized Hamiltonian dynamics. 



^We point out here that the idea of constructing the physical variables entirely in internal terms 
without using any additional gauge conditions is connected with the desire not to distort the 
global properties of the theory and to have all dynamical degrees of freedom under control. 
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non-Abelian constraints by an equivalent set of Abelian ones) and the canonical transfor- 
mation to new variables where a subset of the new momenta is equal to the new Abelian 
constraints. The system of interacting gauge and spinor fields considered in this article 
represent a Hamiltonian system with mixed first and second class constraints. In this 
case the reduction procedure additionally includes the separation of first and second class 
constraints and putting them into the canonical form. 

The paper is organized as follows. In Section II we briefly recall how to obtain the im- 
constrained Hamiltonian system from the initially gauge symmetric one in the framework 
of Dirac constraint theory in order to set the formalism. The Dirac and the Faddeev gauge 
fixing methods as well as the gaugeless method are described. In Section III the gaugeless 
method is exemplified by considering the Yang- Mills system in + 1- dimensions. In Sec- 
tion IV we perform the reduction of the Dirac- Yang-Mills system by explicitly separating 
the first and second class constraints, putting the second class constraints into the canon- 
ical form and Abelianizing the first class constraints. We construct the corresponding 
reduced Hamiltonian system by first eliminating the unphysical gauge degrees of freedom 
and then using the classical scheme of Hamiltonian reduction due to the existence of three 
first integrals of motion. Section V finally gives our conclusions and remarks. 

II. REDUCTION OF CONSTRAINED SYSTEMS WITH FIRST CLASS 

CONSTRAINTS 

The procedure of reduction of phase space of a singular system is a generalization of the 
method of reduction of a system of differential equations possessing a Lie group symmetry. 
The well-known results for this type of reduction in the number of the degrees of freedom 
are embodied in the famous J.Liouville theorem on first integrals in involution. Interest 
in these has revived in connection with the study of Hamiltonian systems with a local 
(gauge) symmetry. Since the work of P. Bergmann and P.A.M. Dirac at the beginning of 
the fifties it has become clear that the role of integrals of motion in a Hamiltonian system 
with gauge symmetry is played by the first class constraints. Although the reduction 
in the number of degrees of freedom due to first class constraints has many features in 
common with the classical case, there are very important differences. In order to explain 
these peculiarities of the reduction procedure and to make the paper self-contained we first 
have to summarize some definitions and to put facts from the Dirac theory of generalized 
Hamiltonian dynamics into the appropriate context. In view of the main purpose of our 
paper, namely to study the finite dimensional system of homogenous Yang-Mills fields, we 
shall discuss the above ideas for a mechanical system, i.e. a system with a finite number 
of degrees of freedom. 

A. The definition of reduced phase space 

Let us consider a system with the 2n - dimensional Euclidean phase space T spanned 
by the canonical coordinates qi and their conjugate momenta pi and endowed with the 
canonical simplectic structure {qi,p^} = Sj. Suppose that the dynamics is constrained to a 
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certain (2n— m) — dimensional submanifold Fc determined by m functionally independent 
constraints 

</'a(p,g) = o, (2.1) 

which we assume to be first class 

{ifiaip, g), ^p{p, q)} = fap^{p, q)^^{p, q) (2.2) 

and complete in the sense that 

Wa{p^ g), Hc{p, q)] = ga^V^ip, q) , (2.3) 

where Hc{p, q) is the canonical Hamiltonian. Due to the presence of these constraints the 
Hamiltonian system admits generalized dynamics described by the extended Poincare- 
Cartan form 

n 

:= E Pidqi - He{p, q)dt (2.4) 

i=l 

with the extended Hamiltonian HE{p,q) that differs from the canonical Hc{p,q) by a 
linear combination of constraints with arbitrary multipliers Ua{t) 

He{p, q) ■= Hcip, q) + Uait)(paip, q) . (2.5) 



^From the condition of completeness ( |2.3|) with He replaced by He it follows that for 
first class constraints the functions Ua{t) can not be fixed in internal terms of the theory. 
This implies that the system possesses a local symmetry and that the coordinates split 
up into two sets, one set whose dynamics is governed in an arbitrary way and another set 



with an uniquely determined behaviour. Recalling the Dirac definition |31| of a physical 
variable as a dynamical variable F with the property 

{F{P, q), Va{p, q)} = dcy{p, q)(py{p, q) , (2.6) 

one can conclude that the first set of coordinates does not affect the physical quantities 
which are defined on some subspace of the constraint surface F^. Indeed, if one considers 
( p.6| ) as a set of m first order linear differential equations for F, then due to the inte- 
grability condition ( |2.2| ) this function can be completely determined by its values in the 
2(n — m) submanifold of its initial conditions [0, 0. This subspace of constraint shell 



represents the reduced phase space F*. This definition of reduced phase space is implicit. 
The main problem is to find the set of 2{n — m) " physical coordinates" Q*, P* that span 
this reduced phase space and pick out the other additional m pairs which have no physical 
significance and represent the pure gauge degrees of freedom. Several approaches to its 
solution are known. Below we shall briefly discuss the corresponding methods of practical 
construction of the physical and the gauge degrees of freedom with and without gauge 
fixing. 
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B. Reduced phase space with the Dirac gauge fixing method 



General principles for imposing gauge fixing constraints onto the canonical variables 
in the Hamiltonian approach were proposed by Dirac in connection with the canonical 
formulation of gravity |3^. According to the Dirac gauge fixing prescription, one starts 
with the introduction of as many new "gauge" constraints 

XaiP,q) = (2.7) 
as there are first class constraints with the requirement 

det\\{xaip,q),MP^Q)}\\ ^ 0. (2.8) 

This allows one to find the unknown Lagrange multipliers Ua{t) from the requirement of 
conservation of the gauge conditions ( p.7|) in time 

Xa = {Xa,Hc} + Y.{Xa,Vf3}Uf3 = (2.9) 

and thus to determine the dynamics of system in a unique manner. A striking result of 
Dirac consists in the observation that such kind of fixation of Lagrange multipliers u{t) 
is equivalent to the following way of proceeding. One can drop both the constraints (|2.1|) 
and the gauge fixing conditions ( |2.7| ) and at the same time achieve the reduction to the 
unconstrained theory by using the Dirac brackets 

{F,G}d := {F,G} - {F,^,}C;X^s',G} , (2.10) 



instead of the Poisson brackets. Here ^ denotes the set of all constraints (^TTI) and (|2.7| ) and 
is the inverse of the Poisson matrix Cap '■= {Ca^^fs} ■ In this method all coordinates 
of the phase space are treated on an equal footing and all information on both initial 
and gauge constraints is absorbed into the Dirac brackets, which describe the effective 
reduction in the number of degrees of freedom from n to n — m 



Y^Ui, Pi, }p.B. = n , ^{gi, Pi, }d.b. 



n — m 



The inclusion of gauge constraints in addition to the initial constraints allows one to 
take the constraint nature of the canonical variables into account by changing the initial 
canonical symplectic structure to a new one defined by the Dirac brackets. The new 
canonical structure, being dependent on the choice of gauge fixing-conditions, is very 
complicated in general and it is not clear how to deal with it, in particular, when we are 
quantizing the theory. However, there is a special case when the Dirac bracket coincides 
with the canonical one and looks like the Poisson bracket for an unconstrained system 
defined on F* 



^Everywhere in the article the dot over the letter denotes the derivative with respect to the 
time variable 
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This representation of the Dirac bracket means that in terms of the conjugate coordinates 
Q*,P* {i = 1, . . . ,n — m) the reduced phase space is parametrized so that constraints 
vanish identically and any function F{p, q) given on the reduced phase space becomes 

F{p,q)\ ^ ^ =F{P*,Q*) ■ 

Thus in the Dirac gauge-fixing method the problem of definition of the "true dynamical 
degrees" of freedom reduces to the problem of a "lucky" choice of the gauge condition. 



C. Reduced phase space with the Faddeev gauge fixing method 

An alternative to the Dirac gauge-fixing procedure has been proposed in the well- 
known paper by L.D. Faddeev [^], devoted to the method of path integral quantization 
of a constrained system. In contrast to the Dirac method, the main idea of the Faddeev 
method is to introduce an explicit parametrization of the reduced phase space. As in the 
Dirac method, one introduces gauge fixing constraints Xa{p,(l) = , but now with the 
additional " Abelian" property 

{Xa(p,g),X/3(p,g)} = 0, (2.12) 

and the requirement ( p.8|) is fulfilled. Now, in accordance with the Abelian character of 



gauge conditions ( p.l2[ ), there exists a canonical transformation to new coordinates 



Qi ^ Qi ■■= Qi {q,p) 

Pi^Pi:=P^{q,p) (2.13) 
such that m of the new P 's coincide with the constraints Xa 

Pa = xAq,p)- (2.14) 



The condition ( p.8| ) allows one to resolve the constraints ( p.lD for the coordinates Qa in 



terms of the {n — m) canonical pairs {Q*,P*), which span the 2{n — m) - dimensional 
surface S determined by the equations 

Pa = 0, 

Qa = QaiQ*,P*) . (2.15) 

After this construction has been carried out, the problem is to prove that the surface S 
coincides with the true reduced phase space F*, independent of the choice of the gauge 
fixing conditions. In other words, it is necessary to find a criterion for gauge conditions to 
be admissible. A radical method to solve this problem is not to use any gauge conditions at 
all. The following subsection will give a brief description of such an alternative gaugeless 
scheme to construct the reduced phase space without using gauge fixing functions. 
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D. The gaugeless method 



If the theory contains only Abehan constraints one can find a parametrization of 
reduced phase space as follows. According to a well-known theorem (see e.g. |^), it is 
always possible to find a canonical transformation to a new set of canonical coordinates 



Pi 



Q^■.-- 

Pi : = 



■ Qi {(i,p) 
Pi{(l,p), 



(2.16) 



such that m of the new momenta, (Pi, . . . , Pm), become equal to the Abelian constraints 

fa 

Pa = Mq,p) ■ (2.17) 

In terms of the new coordinates {Q, P), and {Q*, P*) the canonical equations read 

Q* = {Q*, Hphys} , Q = u{t) , 

P* = {P*,Hp,,y,} , P = 0, (2.18) 

with the physical Hamiltonian 



H,nys{P\Q*) = Hc{P,Q) 



(2.19) 



Hphys depends only on the (n — m) pairs of new gauge invariant canonical coordinates 
(Q*, P*) and the form of the canonical system ( p.l8| ) expresses the explicit separation of 
the phase space into physical and unphysical sectors 



2n { 



Pi 

Qn 
I \PnJ 



2(n — m) 



Q* 
P* 



2m I 



Q 
P 



physical 
variables 

unphysical 
variables 



(2.20) 



The arbitrary functions u{t) enter into that part of the system of equations, which contains 
only the ignorable coordinates and momenta Pa- A straightforward generalization of 
this method to the non-Abelian case is not possible, since the identification of momenta 
with constraints is forbidden due to the non-Abelian character of the constraints. How- 
ever, there exists the possibility of a replacement of the constraints if a by an equivalent 
set of new constraints 



det IIP) I 



ip=0 



^ 



(2.21) 



describing the same surface Tc but forming an Abelian algebra. There are different proofs 
of this statement, based on the resolution of constraints - P, exploiting gauge-fixing 
conditions [^, or using the direct method of constructing the Abelianization matrix as 
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36 



For 



the solution of a certain system of linear first order differential equations 
non-Abelian systems therefore, the construction of the Abelianization matrix and the 
implementation of the above mentioned transformation ( p.l6| ) to the new set of Abelian 
constraint functions $q, completes the reduction of the phase space without using gauge 
fixing functions, solely in internal terms of the theory. 

Before applying the gaugeless method to the construction of the reduced phase space 
of homogeneous Yang-Mills fields in 3+1-dimensional space it seems worth setting forth 
our approach to the same problem in 0+1-dimensional space. 



III. SU{2) YANG-MILLS FIELDS IN 0+1 DIMENSIONS 

In order to explain our main idea how to construct the physical variables we shall start 
with the non-Abelian Christ & Lee model The Lagrangian of this model is 

L := 1 - i V{x^) , (3.1) 

where Xi and i/i are the components of three-dimensional vectors and the covariant deriva- 
tive Dt is defined as 

{Dtx)i := Xi + geijkVjXk . (3.2) 

One can see that this model is nothing else than Yang-Mills theory in + 1 dimensional 
space-time and that is invariant under 5*0(3) gauge transformations. 
Performing the Legendre transformations 

^ I • '^■^) 

dL 

P' = -Q-.=^i+ a^^'^Vj^^k , (3.4) 
one obtains the canonical Hamiltonian 

He = ^PiPi - eijkXjPkVi + V{x^) , (3.5) 

and identifies the three primary constraints Py = as well as the three secondary ones 

$i = etjkXjPk = , (3.6) 

obeying the 5*0(3) algebra 



^In all cases, the proofs use the large freedom in the canonical description of the constrained 
systems. Apart from the ordinary canonical transformations there exist generalized canonical 
transformations i.e., those which preserve the form of all constraints of the theory as well 
as the canonical form of the equations of motion. The Abelianization transformation ( ^.21 ) is 
of course non-canonical, but belongs to this class of generalized canonical transformations. 
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One easily verifies that the secondary constraints are functionally dependent, Xi^., 
We shall now carry out the Abelianization procedure and choose 



,(0) ._ 



:= X3P1 - Xips 



$1 := X2P3 - X3P2 

as the two independent constraints with the algebra 

{$W,$W| = _^$(o)_^$(o) 
X3 X3 



(3.7) 
= 0. 

(3.8) 
(3.9) 



The general iterative scheme of the construction of Abelianization matrix |37] consists of 
two steps for this simple case. Let us at first exclude from the right hand side of eq. 
This can be achieved by performing the transformation 



with the function C obeying the partial differential equation 



X3 X3 

Writing down a particular solution of this equation 

X1X2 



C{x) 



we get the algebra for new constraints 



X3 

Now let us perform the second transformation 



2 I 2 ' 

xi + xi 



^2 ^(1) 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



$(2) .= 

with the function B satisfying the equation 

,(2) X2 



.(1) 



X3 



(3.14) 



(3.15) 



A particular solution of this equation is B{x) = ^ . As result of the two above trans- 
formations, the Abelian constraints equivalent to the initial non-Abelian ones have the 
form 



$1- 



(2) 



X2P3 - X3P2 



X3 



{X3P1 - X1P3) + 



X1X2 

x\ + x\ 



{X2P3 - X3P2) 



(3.16) 
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A. Canonical transformation and reduced Hamiltonian 



We are now ready to perform a canonical transformation to new variables so that two 
new momenta will coincide with the Abelian constraints ( p.l6| ) f\ 

/ 2 I 2 ' P'f' ^2P3 - X3P2 . (3.17) 

\/X2 + Xs 

It is easy to verify that the contact transformation from the Cartesian coordinates to the 
spherical ones 



X2 = r sin sin 6, 9 = arccos 



xi = r cos 9, = V 3^1^ + 2^2^ + a^s^, 

Xi 



\JX\^ + X'^ + Xs^ 



X3 = r cos0sin6', = arctan( — ) , (3.18) 

Vxs/ 

is just the required transformation. Indeed, using the corresponding generating function 
F\x\ pr,pe,P4>] = Pr\/xi^ + + +pe arccos . „ ^\ +p0arctan ( — ) , 

" y/xi^ + X2^ + X3^ ^ VX3/ 

(3.19) 

we get 

sin^ 

PrCos9 — pg , (3.20) 

r 

sm(j)cos9 COS0 

Prsm9sm(p + pe \- P4, — — r , (3-21) 

r r sm 67 

cos0cos^ sin0 

sm cos + — — , (3.22) 

r r sm 9 

and convince ourselves that in terms of these new variables the two independent con- 
straints are indeed pg = and = in accordance with (p.l7|) . It is worth noting 
here that starting with the set of reducible constraints ( p.6| ) and performing the above 
transformation (|3.18|) one obtains the representation 

$1 = -P<p , (3.23) 
$2 = —pe cos (p + P(j> sin cot 9 , (3.24) 
^3 = Pe» sin + p0 cos cot 6' , (3.25) 





OF 


Pi = 


dxi 




dF 


P2 = 


dx2 




dF 


P3 = 


dxs 



^Here we introduce the compact notations for three-dimensional vectors x,p and multiply 
the constraint $2'^^ by the factor \/x2^ + xs'^ to deal with constraints of one and the same 
dimension. This multiplication conserves the Abelian character of the constraints, since 

{^f\Vx2'+X3'^} = 0. 
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adapted to the Abehanization. The corresponding Abehanization matrix for the reducible 
set of constraints is 



^ / —6/2 sin (f) — cos , di sin , di cos , \ 

D := - ((i2 COS0 — (is sin 0) cot 6* , — (is — rfi cos cot 6^ , (i2 + (ii sin cot , , (3.26) 
\ cot 6* , sin , cos , / 

with arbitrary d and d := di cot 6 + (i2sin0 + (is cos 0. This example demonstrates two 
important features of the Abehanization procedure: it is not necessary to work with an 
irreducible set of constraints, because the Abehanization procedure leads automatically 
to an irreducible set of constraints, ii) in certain special coordinates the problem of the 
solution of differential equations reduces to the solution of a simple algebraic problem. In 
terms of the new canonical variables the canonical Hamiltonian (|3.5|) reads 




He = -Pr + 7^ {Pe + ^ - Pm - Poye + V{r) , (3.27) 



with the physical momentum Pr = —, — ^2^^ , and 

y<i>-=yi + y2 sin + 1/3 cos cot 61 , 

1/0 := I/2COS0 - ?/3sin0. (3.28) 

As a result, all the unphysical variables are separated from the physical r and and 

their dynamics is governed by the physical Hamiltonian obtained from the canonical one 
by putting p^ and pe in ( p. 271) equal to zero 

Hphys = \pI + V{r) . (3.29) 



IV. SPATIALLY HOMOGENEOUS SU{2) DIRAC-YANG-MILLS FIELDS IN 

3+1 DIMENSIONS 

A. Canonical formulation of the model 

The dynamics of SU{2) Yang-Mills gauge fields A^(a;) minimally coupled to the 
isospinor fields \E'q,(x) Q in four-dimensional Minkowski space-time is defined by the La- 
grange density 

^ = ^Y.-M. + ^Matter + ■ (4.1) 

The first term is the kinetic term of the non-Abelian fields 



The matter isospinor variables are treated classically as a collection of four Grassmann 
quantities. Detailed notations are collected in the Appendix. 
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Cy.-m. = ^ tr (F^.F^'^) , (4.2) 
the second term corresponds to the matter part 

C-Matter = ^ [^a^ft^d^^c. " S'^^aT/.^a] " ^ ^ a'^ a , (4.3) 

and the last term describes the interaction between the gauge and the matter fields 

Cr = g ^^',7^(r'^)„^vl>^ A^^ , (4.4) 

with the Pauli matrices Ta, a = 1, 2, 3. 

After the supposition of the spatial homogeneity of the fields, (^]l|) reduces to a finite 
dimensional model described by the Lagrangian 

L = ^{DtA)ai{DtA)ai + ^ (^a70^a " ^a70^a) " - QPa Ya + QJiaAai - V {A) , 

(4.5) 

where the nine spatial components A^ are written in the form of a 3 x 3 matrix Aai-, the 
time component of the gauge potential is identified with Ya := Aq and Dt denotes the 
covariant derivative 

{DtA)ai '■= Aai ~ Q^abcYbAci ■ 

The part of the Lagrangian density corresponding to the selfinteraction of the gauge fields 
is gathered in the "potential" V{A) 

V{A) := ^ [tr'(AA^) - ti{AA^f] , (4.6) 

while their interactions with the matter fields are via the isospinor currents 

1 , 



jiam := ^^ali{ra),p^p . (4.7) 

After Legendre transformation one obtains the canonical Hamiltonian 

He = ^EaiEai + m^^<i/^ - g {eabcAdEbi - Pa) Ya - gjiaAai + V{A) , (4.8) 

defined on the phase space endowed with the canonical symplectic structure (see Ap- 
pendix) and spanned by the bosonic and fermionic canonical variables {Ya, Pva), {Aai, Eai) 
and {^!a,P^^), (^a,P^^), where 

r)T 

Py.:=-^ = 0, (4.9) 

oYa 

dL 

Eai ■■= -TTT- = ^<^i - aeabcYbAci , (4.10) 

dAai 

P^^:=L^ = -'-^.7o, (4.11) 
P^. ■■= tJ-^ = -^7o^a . (4.12) 
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According to the definition of tlie canonical momenta ( [4.9|) , ( |4.11| ) and ( |4.12|) tlie pliase 
space is restricted by tlie three primary bosonic constraints 



P-(^ = , (4.13) 
and the sixteen Grassmann constraints 

:= P*. + '-^.10 , := P^^ + . (4.14) 
Thus the evolution of the system is governed by the total Hamiltonian 

Ht := He + uUt) + ulit) + ulit) . (4.15) 

The conservation of bosonic constraints (|4.13| ) in time entails the following further con- 
dition on canonical variables 



^a■.= eabcAc^Eu- PaM=^ , (4.16) 



which is the non-Abelian Gauss law. In contrast, the maintenance of Grassmann con- 
straints and in time allows to determine the Lagrange multipliers u\{t) and M^(t) 
in the expression ( [4. IS] ) for the total Hamiltonian. Taking into account the Poisson brack- 
ets of constraints 

$j} = eijk^j + e.jkpk [^] , (4.17) 
T^} = -^^70 (ra )^„ , (4.18) 
{$,,T2}= 7o(ra),^vl/^, (4.19) 
{T^, Tj} = -i5^p 70 , (4.20) 

one can convince oneself that no new constraints emerge and hence that ternary con- 
straints are absent in the theory, $ = . 

Constraint Shell 

To implement the reduction procedure without using gauge fixing conditions we have 
to put the constraints into the canonical form discussed in the next paragraph. 



B. Putting the constraints into the canonical form 

1. Separation of first and second class constraints 



= 16. Hence among the constraints there are six first 

Ca = 



The set of the 22 constraints := (Py,$,T) represent a mixed system of first 
and second class constraints. The Poisson matrix M^g := {C_4, Cg} is degenerate on 

constraint shell, rank 1 1 1 1 

class ones. 

In order to perform the reduction procedure let us start with the separation of the first 
and second class constraints. The primary constraints Py "commute" with all the other 
constraints and thus we should deal only with constraints := ($, T). The separation of 
constraints is achieved by a transformation to an equivalent set of constraints := ($, T) 
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C'a--=D'^bC's. (4.21) 

so that the first class constraints $ form the ideal of the algebra 

{l>,T} = 0, {<|,|.}|^^,= 0, (4.22) 

and the pairs of second class constraint satisfy the canonical algebra 

{Ti^Tl} = - . (4.23) 

In order to transform the algebra of constraints to the canonical form let us at first perform 
the equivalence transformation 

$1 := + Tl'- (r,)„^ + (r.)^„ , (4-24) 
on the bosonic constraints and the equivalence transformation 

t^=-zT;,7o, ^1-.= rl, (4.25) 
on the Grassmann constraints. The Poisson brackets of the new constraints 

{$1, <l>a = eatcK , (4.26) 



{rlK}= ;^Tj(rJ,„, (4.27) 



{TL<^>l} = -^(Ta)„,Tj, (4.28) 
{flfl} = -t6^p , (4.29) 

show the separation of the first class constraints on the surface T = defined by the second 
class constraints. To achieve this separation on the whole phase space it is necessary to 
apply the additional transformation 

^a:=K-mra)^,fl. (4.30) 
One can verify that the first class constraints form the ideal of the total set of constraints 

$b} = e,bA, (4.31) 
{t;,,l>4 = 0, (4.32) 
{fl^a} = (4.33) 

and the second class constraints obey the canonical algebra ( |4.29| ). The explicit form of 
the resulting set of constraints is 

Ti :=-^P*„7o + ^vI/ = 0, (4.34) 

^':=^*. + ^ 70^ = 0, (4.35) 

$a := eabcAb^E^^ + i^'r„7o^ - ^P^Ta-foPq, + ^ [P^Ta"^ + ^T^P^) = . (4.36) 
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In order to implement the reduction due to the second class constraints ([4.34]) and ( |4.35|) 
let us introduce the new canonical variables and (11^^, n^^) via 



1 / — 



n 



7o 



In terms of the new variables the constraints read 



(4.37) 
(4.38) 



a 
-2 



= ^abc^biEci — -11^^ TaQ^^ = 



(4.39) 

(4.40) 
(4.41) 



2. Canonical transformation to adapted coordinates 



The example of the Christ and Lee model in Section |TI| shows that the realization of 
constraints by Abelianization is immediate if one performs a canonical transformation to 
a new set of variables containing the gauge invariant ones as a subset. Hence in order to 
simplify the Abelianization of constraints let us single out the part of the gauge poten- 
tials Aai, which is invariant under gauge transformations. Because under a homogeneous 
gauge transformation the gauge potentials transforms homogeneously one can achieve the 
separation of gauge degrees of freedom by the following simple transformation 



Aai{Q,Q*)=Oak (q) Q- 



ki 1 



(4.42) 



where O is an orthogonal matrix, O G SO (3), and Q* is a positive definite symmetric 
matrix. This transformation induces a point canonical transformation linear in the new 



canonical momenta. The new canonical momenta (P*^ 
generating function 



Pi) can be obtained using the 



F4 [E- Q,Q*)=J2 Ea^Aa^ (Q, Q*) = tr iO{Q)Q*E'^ 



(4.43) 



as 



p* 



ik 



Wlk 2 



tr 



- (OE^ + EO^) . 
2 V Afe 



(4.44) 
(4.45) 



In order to express the Hamiltonian and the Gauss law constraints in terms of these 
new canonical pairs let us write the field strengths Eai in the form 



Eai — Oak (q) Lki (P, P*; Q, Q* 



(4.46) 
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with a 3 X 3 matrix L^i to be determined. One can immediately see that the symmetric 
part of the matrix L is equal to the new momenta P* 



P*ik — o (-^'fc + ^ki) 



(4.47) 



and a straightforward calculation shows that its antisymmetric part is 



■7: {Lik — L/^i) — ej/fc (7 ^)ls ^) . Pj — (-msn {P* Q 



with 



and 



0{Q) 



(4.48) 



(4.49) 



7ik ■■= Q*ik - Sik tr(Q*) . 



(4.50) 



Thus the final expressions for field strengths Eai in terms of the new canonical variables 
are 



Eai = 0{Q)ak 



(4.51) 



3. Abelianization of first class constraints 

The formulation of the theory in terms of the new variables is adapted to the procedure 
of Abelianization. Using the representations ( |4.42| ) and ([4.51| ) one can easily convince 
oneself that the variables Q* and P* make no contribution to the secondary constraints 
( [4.41| ) and Q , P enter well-separated from the physical matter variables 

:= OUQ)n~^,Pj - ^n^„ ir,)^pQl, = . (4.52) 

In order to deal with the Abelianization it is useful to perform the following canonical 
transformation on the Grassmann variables 



=■■ n.Upa (Q) , (4.53) 
=■ (g) Ql, . (4.54) 

with the unitary matrix U in the two dimensional representation of 5*0(3) chosen such 
that 

Oab = ^tr (u+TaUn) . (4.55) 
As a result, the Gauss law constraints ([4.52|) take the form 
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$1 := ^l-\P, - Vl^ (rj„^ Q*, = . (4.56) 

Hence it is clear that the matrix VL~^ is just the matrix of Abehanization D in ( |2.21| ). 
Hence, after performing the Dirac transformation with the matrix D := Vt{Q) on the 
constraints the equivalent set of Abelian constraints is 

Pa-^ asQs = , (4.57) 

with 

e.:=^n.(r,)„,Q^^. (4.58) 



C. Reduction due to the Gauss law and the second class constraints 

In the previous section, in accordance with the general scheme of reduction formulated 
in subsection ( |11D| ), the new set of constraints in canonical form have been obtained and 
the adapted canonical pairs been chosen for the explicit implementation of the Gauss laws 
( [4.16| ) and the second class constraints ( [4.29|) . After having rewritten the model in this 
form, the construction of the unconstrained Hamiltonian system is straightforward. In all 
expressions we can simply put P = QQ and H^^ = Q^^ = 0. In particular, in terms of 
the "physical" electric field strength £ai 



=■■ OakiQ) £M*,P*) , (4.59) 

p=ne 



the physical unconstrained Hamiltonian 

Hphys '■= Hc{P,Q) 

may be written as 



constraint shell 



^5,-/-*'- = ltT{S') + ^-^[tT\Qr-tTiQr] -^7tr(j*g*)+^m(n.7oQk) • (4-60) 
where j* is the isospin current in terms of the new Grassmann variables 

jI:=^n.7.7o(ra)„^Q^,. (4.61) 



With the aid of the identity det7 eisk {l~^)si = ^aiblialkb and representation ( |4.51|) 
for the field strengths, we find the explicit form for the "physical" electric field strengths 
in terms of P* and Q* 

Sk.{Q\ P*) = P:, + -^eakiT (^MrJi) , (4.62) 
det 7 

where Ai denotes the isospin angular momentum tensor 
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Mmn ■= emsnJs ■ (4.63) 

Here 

J,:=e, + T,. (4.64) 

is the sum of the gauge field isospin vector := {Q*P*)mn matter field 9^ 

defined in ( [4.58| ). With ( |4.62| ) the unconstrained Dirac- Yang-Mills Hamiltonian reads 



2 



H^^I-^- = §tr(P*)^ + ^^tr i^M^r + J [tr\Qr - tr(Q-^ 

-gtTifQ*)+zm(vi^oQl) . (4.65) 



In order to achieve a more transparent form for the reduced Dirac- Yang-Mills system 
( [4.65| ) one can perform a canonical transformation expressing the physical coordinates 



Q* and P* in terms of new variables adapted for the analysis of the rigid symmetry 
possessed by the reduced Hamiltonian system ([4.65|). It is convenient to decompose the 



nondegenerate symmetric matrix Q* in the following way: 

Q* = n^{^,e,(P)vn{ip,e,(i)), (4.66) 

with the 50(3) matrix TZ parametrized by the three Euler angles Xi '■= (V'? 0); (see Ap- 
pendix) and with the diagonal matrix T> := diag (xi, X2, x^) . The corresponding canonical 
conjugate coordinates {Pip,Pe,P<f>,Pi) can be found by using the generating function 

F [xi, iP, 9, 0; P*] := tr (Q*P*) = tr (7^^(x)/^(x)7^(x)P j (4.67) 

as 

OF 

Pi=Q^.= (P*7^'^a^7^) , 

BF ( (YR \ 
Vx. = ^ = tr y-^n \P'Q' - j , (4.68) 

where a, are the diagonal members of the orthogonal basis for symmetric matrices = 
(oj, ai) 2 = 1,2,3 given explicitly in the Appendix. The original physical momenta P*^ 
can then be expressed in terms of the new canonical variables as 



\s=l s=l / 



(4.69) 



with 



Vi := — ^ — , {cyclic permutation i ^ j k) (4.70) 

Xj Xjf 



and the 5*0(3) left-invariant Killing vectors 
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sin ijj 

:= -^—^ Ps + cost/' pe - sin?/' cot 61 p^ , (4.71) 
sm 

cos ip 

6 := P4> + sin^ + cos?/' cot 9 p^ , (4.72) 

sm 

6:=^^. (4.73) 
Representing the physical electric field strengths Sai in the alternative form 

= P\k + {lJsl),k Js , (4.74) 

det 7 

with the SO (3) generators Js given in explicit form in the Appendix, we finally get the 
following physical Hamiltonian defined on the unconstrained phase space 

P'^^^ 2^^^ + 4^'^ + 4 ^ U,+a;J ^2 ^ * ^ 

s=l s=l cyclic \ J / i<j 

-gjZj^s + imivi^l.Ql}) . (4.75) 

s=l 



Note that for the pure Yang-Mills system ( 4.75 ) reduces to 



s=l ^ cyclic yXj — XI j i<j 

This completes our reduction of the spatially homogeneous constraint Dirac- Yang-Mills 
system to the equivalent unconstrained system describing the dynamics of the physical 
dynamical degrees of freedom. However, apart from this reduction due to the underlying 
gauge symmetry, there is the possibility to realize another type of reduction connected 
with the rigid symmetry admitted by the unconstrained system (^4.65 ). For simplicity, 
the discussion in the next section will be restricted to the pure Yang-Mills system and 
we shall show how to further reduce the obtained 12- dimensional system ( ^.76| ) to an 
8-dimensional one in general and, for a special case to a 6-dimensional one using the 
corresponding first integrals. 



D. Further reduction using first integrals 



The reduced Yang-Mills theory ([4.76|) has a rigid symmetry connected with the exis- 
tence of the first integrals 



^ijk 



ak ■ 



(4.77) 



For the subsequent reduction in the number of degrees of freedom we shall use the integrals 
of motion ( 4.77 ). One can verify that in terms of the new variables they read Jj = TZik^k ■ 
In contrast to the reduction due to first class constraints the values of first integrals are 
arbitrary and depend on the initial conditions. In this case reduction means to consider 
the subspaces of phase space which are the levels of fixed values for these first integrals 
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(4.78) 



and the subsequent construction of the quotient space with respect to the rigid symmetry 
group. Therefore, in contrast to the reduction that we have done before, the first integrals 
in general are a mixed system of first and second class constraints. By "in general" we 
mean that not all constants of motion Cj are zero. In this case the rank of the Poisson 
matrix is rank | | ~ ^ ' which means that there is one first class constraint 

and one pair of second class constraints []. The problem is now to separate the algebra of 
constraints and to find the equivalent set of constraints \E'i = , so that the first class 
constraint \E'i forms the center of the algebra 

{^l,^^} = (4.79) 

and the pair of second class constraints obey the canonical algebra 

{^2,^3} = 1. (4.80) 

After having passed to new variables in the last section in order to isolate the gauge degrees 
of freedom from the physical ones, we shall now perform another canonical transformation 
from the physical variables to new physical variables so that one of the new momenta 
coincides with the first class constraint \E'i and another pair of new canonical variables 
coincides with the pair of second class constraints \l/2 and \l/3. In terms of these new 
canonical variables the reduced system is obtained by reducing the Hamiltonian ( |4.65| ) to 



the integral surface ( [4.78|) . As result the new Hamiltonian will depend on 4 canonical pairs 



and one parameter which reflects the existence of the integrals of motion. To demonstrate 
this let us choose the integral constants as q = (0, 0, c) without loss of generality. One 
can then write down the needed new set of constraints, describing the surface ( |4.78D , in 
the form 



:= If + + l^^-c' = 

'h 



^2 ■■= arctan ( :^ ) = , (4.81) 



^3 := /i = . 

We are now ready to perform the transformation to special canonical variables so that 
the pair of second class constraints is equal to the one pair of the canonical variables and 
one equal to the new momentum Q 

Ho := ^1 , Hi := ^3 , X, := ^2 • (4.82) 

and complete the set of canonical coordinates by the following pair 



^In the exceptional case when Cj = we can consider the three integrals as first class constraints, 
and this circumstance leads to a further reduction of our system. 

^ This type of variables are well-known from rigid body theory as Depri 
variables used in celestial mechanics. 



41] or Andoyer m.2c\ 
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Xo := arctan 



.6 



(4.83) 



The canonical conjugate coordinate Xq can be determined with the help of the generating 
function 



- — Jul sin^ a - nf - ni + 2nin2 cos 

sma ^ 



a . 



(4.84) 



Due to the symmetry the Xq is a cyclic coordinate and the reduced Hamiltonian depends 
only on the canonical pair 



:= arctan 



.6 



— 



Hence, using the first integrals ( |4.77] ), the pure Yang-Mills Hamiltonian ( ^.76| ) can be 
further reduced to 



ttY.-M. * 

phys 



1 ^ 1 
^ s=l ^ 



2 - sm ^— -2 - COS l\} 



(/-y>2 ,-y>2 



(xi - a;i) 



X2 ^2 \ ^ 
3 -^iJ J 



2 



(4.85) 



where in accordance with the general scheme of reduction there arises the additional 
so-called reduced potential term 



X2 I ,y»2 /y>2 I /vj2 

sm"?^^^f l^ + cos^V' 



(a^i - a;!) 



(xi - xf ) 



2^2 



(4.86) 



It is interesting to point out the difference between the reduced Yang-Mills Hamiltonian 
( [4.76| ) and the corresponding one in the recent work by B. Dahmen and B. Raabe. In 
contrast with their representation for the gauge potentials, in which the gauge degrees of 
freedom are mixed with the rigid rotational cyclic coordinates, we have started with the 
explicit separation of all physical degrees, including the rotational ones. And only after 
the reduction in the number of degrees of freedom due to the rigid symmetry the obtained 
Hamiltonian ( [4. 85] ) coincides with the one obtained in the work by B. Dahmen and B. 



Raabe |2J] for pure Yang-Mills mechanics. 



V. CONCLUDING REMARKS 



As mentioned in the introduction our investigation has pursued two goals. One is 
pure theoretical interest. Due to the homogeneity condition SU(2) Dirac- Yang-Mills field 
theory has greatly simplified to a finite dimensional mechanical system, for which one can 
describe the equivalent unconstrained system in an explicit way. However, apart from this 
reason, there is also an interesting application of this model. It has been known for a long 
time, that, if one considers the Euclidean QCD effective action as a function of the non- 
abelian electric and magnetic fields E and 5, one finds that there are field configurations. 
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corresponding to nonvanishing E and B fields, for wliicfi tfie value of the effective action is 



lower than that for E = and B = |^0[. This observation indicates a drastic difference 
between the true ground state of QCD and the corresponding perturbative vacuum and 
constitutes the basis of all models of condensates. One of the main reasons to study the 
dynamics of spatially constant Yang-Mills fields, is the faith that the corresponding zero 
momentum quantum operators are very important for the description of the QCD ground 
state due to the presence of the IR singularity. There are many attempts to exploit 
the homogeneity approximation for gluon fields with the aim to shed light on the vacuum 
structure of QCD. We also adhere to this position and our task in this note was to prepare 
the classical description of Yang-Mills mechanics in a form that we are going to exploit 



for the description of squeezed vacuum 43 
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APPENDIX A: NOTATIONS AND SOME FORMULAE 



1. Definition of configuration variables 

SU{2) Dirac- Yang-Mills theory considered in this paper includes as dynamical vari- 
ables the set of spin-1 gauge fields := ^4^x^/2, a = 1, 2, 3 in the adjoint representation 
of SU(2), with the corresponding field strengths 

■■= r72, (Al) 
F;, := d,Al - d,A1 + ge'^^'^AlAl , (A2) 

and the matter spinor (Dirac conjugate spinor) field variables in the fundamental 

representation of SU(2) with values \E'a := (\E'^, . . . , obeying the Grassmann algebra 

^>'^J>^, + = . (A3) 



2. Hamiltonian structures 



Generalized Poisson brackets for functions on a phase space spanned by both even and 
odd coordinates Za '■= {{Y, Py), {Ai, Ei); (\I/q, P^^)) are defined as 
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{F{Z),G{Z)} ■.= Y.F^^AB^G . 

The nonvanishing components of the canonical symplectic form ujab '■- 
expUcitly 

{y„,P^} = (5^ {A,,, E'^}^ SIS', 

for bosonic degrees of freedom 



(A4) 

{Za,Zb} read 
(A5) 



(A6) 



for fermionic degrees of freedom. 



3. The Euler parametrization for S0(3) group 

The conventional representation of SO {3) group elements in terms of Euler angles 

has been used in main text with the following matrix realization for the generators Jj 
obeying the 5*0(3) algebra [J,, Jj] = eijk Ju 



(A8) 





/O \ 




/O 






/O 


-1 




Jl = 


0-1 











1 









lo 1 y 




i,l 














4. Basis for symmetric matrices 

We use the orthogonal basis = a*) for symmetric matrices. They read explic- 



itly 



1 





a2 




1 




Q!3 





1 



(A9) 





/O 0\ 




/O 1 \ 


/O 1 0\ 




1 


, — 


0, = 


1 




lo 1 oj 




li o; 


i^O 0/ 



They obey the following orthonormality relations: 
tr (cfjCfj) = Sij , tr (cCjCCj) = 2Sij , 



tr (aiaj) — . 



(AlO) 



(All) 
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